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ABSTRACT 


Selective Reenlistment Bonuses (SRBs) are offered to 
improve retention in designated military occupational 
Specialties (MOSs) for specified years-of-service intervals 
(zones). The amount of the bonus is set by assigning an "SRB 
Multiplier" for each MOS and zone combination (cell). 

Determination of multipliers is modeled as a nonlinear 
knapsack problem which is then linearized to a generalized 
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provides lower bounds and feasible solutions. The best 
feasible solution is improved uSing a greedy heuristic to 
apportion unexpended funds. 

A FORTRAN 7/7 computer program implements the procedure. 
Data for FY86 yields a O-1l integer program with 4895 binary 
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I. INTRODUCTION 


The United States Marine Corps (USMC) uses Selective 
Reenlistment Bonuses (SRBs) to increase retention in desig- 
nated military occupational specialties (MOSs) and years-of- 
service intervals (zones). Those MOSs and zones which will 
receive bonuses as well as the dollar amounts of the bonuses 
are determined by the assignment of "SRB Multipliers.” In 
this study, the problem of assigning multipliers given a 
limited budget is modeled and solved as a mathematical 
program. In this chapter, the SRB program is described, 
terms are explained and the problem itself is defined. In 
addition, the solution approach is briefly discussed and the 


structure of this thesis is outlined. 


A. PROBLEM STATEMENT 

The Marine Corps SRB program was developed as a reen- 
listment incentive to increase retention in designated MOSs. 
In effect, an SRB is a sum of money offered to individuals 
with specific skills and years in service to encourage them 
to reenlist and thereby selectively increase manning levels. 
The bonus is applied at reenlistment points that fall 
between 21 months and 14 years of active duty service (ADS). 
This period is divided into three reenlistment zones. As 
defined in the applicable Department of Defense (DoD) direc- 
tive [Ref. 1], the period between 21 months and 6 years of 
ADS is designated as Zone A, the period between 6 and 10 
years as Zone B and the period between 10 and 14 years as 
Zone C. 

SRBs are not offered to everyone who reenlists, rather 
they are directed at selected MOS and zone combinations, 
called "cells" in this study. For each cell, there exists a 


reenlistment target which is defined as the desired number 


of reenlistments in that cell. In many cells, enough Marines 
will reenlist without any added inducement so that the reen- 
listment targets will be automatically met or exceeded. In 
others, however, the absence of a bonus will result in 
insufficient reenlistments. Through the offer of bonuses to 
reenlist, it is possible to reach or approach the targets. 
To permit additional flexibility in awarding bonuses, the 
amount of an SRB can vary from a minimum of 0 to a maximum 
omer s LO. 000). From experience, it is known that larger 
bonuses result in greater retention rates although the rate 
at which retention increases with the size of the bonus 
differs greatly between cells. The total amount of bonus 
money that can be awarded per fiscal year is limited by the 
SRB program's annual budget. Since an individual receives 
(Lat of the bonus as a lump sum payment with the remainder 
apportioned in equal annual payments over the term of the 
reenlistment, the total SRB budget must consist of two 
parts. The first is money owed from bonuses from previous 
years and its apportionment is non-discretionary. The rele- 
vant budget here is the remainder which can be allocated for 
current year bonuses. 

Once the decision is made to consider the use of a non- 
zero SRB to improve manning in a particular cell, the actual 
monetary size of the bonus is based on increments of the 
eligible member's monthly base pay. The individual receives 
an amount equal to his monthly base pay multiplied by that 
cell's SRB multiplier and the number of years of additional 
obligated service, subject to a maximum of $16,000. The DoD 
directive stipulates that the multipliers may not exceed 6 
and may be fractional amounts. The Marine Corps’ implementa- 
tion of this policy is somewhat more restrictive. 1Gis 
chooses not to uSe fractional multipliers and does not 
permit a multiplier greater than 5 for reenlistments in Zone 
ae esemembe ane 5 in)Zone C. Assigning a multiplier of 0 


1s equivalent to no bonus. 


Presently, SRB multipliers are assigned manually in an 
lterative procedure. Multipliers are xP eg cells 
based on the desired number of reenlistments for those 
cells. Once a particular combination of SRB multipliers has 
been selected, a microcomputer based spreadsheet package is 
used to determine if the proposed combination is feasible, 
1.e., within the allowed budget. If it is not, multipliers 
are changed and the calculation is repeated: 

At this time, there is no objective function with which 
alternative sets of proposed multipliers can be evaluated or 
compared. Instead, the officer responsible for assigning 
multipliers Penies exclusively on his experience, judgement 
and Various <rudes of «thumb tomyucdee athe desirabi li liye 2 
particu lagasSe of multipiiersminvenmie alsensible wobyecemee 
function were defined, it is doubtful that an optimal solu- 
tion could be determined uSing a manual procedure. Of the 
approximately 1000 cells, on average, 450 are assigned non- 
zero SRB multipliers during any given period. It. 1S @gemy 
likely that the dimensionality of the problem would quickly 
overwhelm any ability to manually arrive at an optimal 
allocation. 

This study proposes an alternative procedure. It is 
automated due to the large amount of data. An objective 
function is defined and justified, thereby permitting the 
comparison of different sets of multipliers. Further, an 
optimization method built around this objective function and 
the budget constraint is described and tested. Using this, a 
nearly optimal set of SRB multipliers can be determined. 
This will result in a more effective allocation of the SRB 


Program budget. 


B. PROBEEMPSCOrE 
The goal of the SRB program is to reduce expected short- 
falls in the number of reenlistments in particular cells by 


imcreaSing the retention rate in those cells through the 


offer of a reenlistment bonus. The model that is developed 
in this study tries to achieve that goal by determining a 
set of SRB multipliers which, by a selected definition , can 
be considered optimal. This 1s a single period model; no 
attempt is made to forecast sets of multipliers for subseq- 
uent years based on expected retention resulting from the 
current year's set of multipliers. It allocates that part 
of the budget not obligated for annual payments for bonuses 
from previous years and does not consider the effect of the 
257, of the bonus that is apportioned in future years. 
Furthermore, the model, through its use of expected values, 
is essentially deterministic. Other assumptions regarding 
the model are explicitly addressed in the succeeding 


chapter. 


C. MODEL AND SOLUTION APPROACH 

In the mathematical formulation that is proposed to 
model the assignment of SRB multipliers, the objective func- 
tion to be minimized is the sum of individual cell penal- 
ties. These penalties are nonlinear increasing functions in 
the deviation from the reenlistment target for the cell. 
Inclusion of the budget constraint, which is also nonlinear, 
results in a nonlinear knapsack problem. The nonlinearity, 
coupled with other complicating factors, prevents’ the 
problem from being efficiently solved using common knapsack 
algorithms. For descriptive and computational purposes, the 
problem is reformulated as a special generalized assignment 
problem. 

In this generalized assignment formulation, cells repre- 
Sent the tasks that must be assigned. Unlike conventional 
generalized assignment problems, there is but a single agent 
who must accomplish all the tasks but who has alternative 
processes for each. These alternative processes correspond 


to the various multiplier levels permitted for that cell. 


Because of the special structure of this_' problem, 
Lagrangian relaxation is particularly well suited as a solu- 
tion technique. Lagrangian relaxation has been successfully 
used in scheduling (e.g., [Reta location 
(e.g., [Ref. 3]) and set-covering problems (e.g., [Ref. 4]) 
as well as generalized assignment problems [Ref. 9]. Of 
interest is that previous applications to generalized 
assignment problems have shared the characteristic of 
multiple tasks and multiple agents; a review of the litera- 
ture has not revealed an application with a formulation 
Similar to the one developed in this study. 

An important advantage to the Lagrangian relaxation 
approach is that, in the process of establishing a lower 
bound on the optimal solution, feasible solutions are also 
obtained. A heuristic method is developed which improves 
the best feasible solution uncovered in the bounding proce- 
dure. The heuristically improved set of multipliers is 
accepted as the final solution provided its objective value 


is sufficiently close to the lower bound. 


Do THESTSSOUTEENE 

This thesis develops and presents a method for deter- 
mining SRB multipliers in the USMC. In Chapter II, the 
problem is formulated first as a nonlinear knapsack problem. 
Then, to facilitate the description of the solution approach 
and for computational purposes, it is reformulated as a 
Special generalized assignment problem. The solution meth- 
odology is presented in Chapter III. The approach uses the 
technique of Lagrangian relaxation combined with a heuristic 
procedure to provide a final solution. Details on the 
implementation of this process and computational results are 
presented in Chapter IV. Conclusions and recommendations 
are contained in Chapter V. Lastly, a listing of the source 
code, uSer instructions and samples from input and output 


files are included as appendices. 
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Pine rORMeEATION AND DESCRIPTION 


In this chapter, the process of determining multipliers 
for all cells is formulated first as a simple nonlinear 
integer program. Coefficients are defined and described and 
the rationale behind the selection of the objective function 
is explained. For explanatory and computational purposes, 
the problem is then converted to a special case of the 


generalized assignment problem. 


MODEL. DEVELOPMENT 
The following formulation is developed to model the 


assignment of SRB multipliers: 


EnGaee's : 
ae ee cells 
Data. 
B budget 
Bumet tons: 
Cee ) penalty for setting multiplier at x; 
iemeree 1.1. 
r; (x; ) ecto sclbune multiplier at x. 
in cell 1 
na maximum permissible multiplier for 
Sete. = DS oeeeetom Zone A. 8, C 
respectively 
Decision Variables: 
Xi Miter plwer for cell 1 
Formulation: 
min YOHES (P1) 


Se Dri (sy) =a 
{ 


i 


= i= 1 een 


Sj LINE CIES 


In defining the objective function it is assumed that 
each cell’s effect on the overall” objective |funeerion sis 
independent of all other cells. The objective function, 
then, is just the sum of individual cell penalties where a 
c.(x.) should be 0 Z£ the desired number (of Maj yos ame 
expected to reenlist given an SKB multiplier =els eee 
Should rise as the expected deviation from the reenlistment 
target increases. An individual term in the budget 
constraint represents the expected reSource requirements 
associated with selecting a particular multiplier value for 
that cell: For this formulation the following functional 


formsSeare Seleecrea. 


ea (xs) = Graig ala) (eqn 2) 
Ch@=oRs Gene ale bs eGR GoyyE.-D;) 1°} 


Be ch = (R; (x; JE; ) {min Neva Eee) 12000) } (eqn 2 ay 


where 
a normalized cost of training a Marine in cell i 
a number of Marines in cell i 
W; exogenous cell weighting factor 
De reenlistment target for cell i 
R,(x,) reenlistment Tesponse prate  inec oleae a vhem 

multiplier x. Sasgethened 

Ee number of Marines in cell i eligible for reenlistment 
Ee average monthly base pay of a Marine in cell i 
Si duration of reenlistment of a Marine in cell i 
Q relative weight of overages to shortages with 


respect to the reenlistment target 


(a-b)” positive part of (a-b) 


WW 


The function c;(x;) represents a weighted squared devia- 
tion from the target number of reenlistments that would be 
Spe eeedeimncell yi eit mWitiplier x; were offered. It is the 


product of two sets of terms. The first term ie. ow repre- 


al 
Sents a "weight" assigned to each cell. A weighting scheme 
is required as deviations in some cells are considered more 


critical than in others. The second set of terms 
(( (Ds -R; (x; )E,)*]° - QL (R; (x; )E,-D,)*17) 


is the squared deviation from the desired number of reen- 
listments in cell i but where Q may be used to weight the 
penalty incurred by overages (the second difference term) 
differently from the base penalty incurred by shortages (the 
Mieet ditterence term). ThewQ, the "over/under’ “factor, has 
been set to l in this study but values somewhat less than l 
might be appropriate. 

Squaring the deviation has the effect of penalizing 
large deviations much more severely than smaller ones. This 
is appealing because it evens out shortages and overages 
among cells. Most would agree that it is less detrimental 
to the force as a whole to "spread the shortages around” 
than to have a potentially debilitating shortfall in a few 
cells. Spreading overages around is desireable for logis- 
tical reasons. 

The deviation terms are calculated by squaring the 


difference between the desired number of reenlistments D. 


and the expected number of reenlistments Re Coane The 
function R,;(x;) is discrete and defined only at integer 
Vale ote son the interval irr eemopccitacally, whensthe 


individuals in cell i are offered an SRB multiplier x a 


5 ae, 
fraction R;(x,) of the E,; individuals eligible to reenlist 
will Choose to do so. The values of R,(x;) have been esti- 
mated from responses to previous bonus offers; more sophis- 


ticated methods for determining these functional values, for 
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example, by including the effect of past and present 
economic conditions, are currently being studied [Ref. 5]. 

As stated previously, otherwise identical deviations in 
different cells are not viewed equally. Generally it costs 
the Marine Corps considerably more to replace or train an 
individual in a highly technical MOS than in a less tech- 
nical MOS. Further, the actual number of Marines in a cell 
influences the acceptability of shortages. For example, a 
shortfall of 5 in a cGell containing 10007 could {be mene 
easily tolerated than that same shortfall in a cell 
COnEadinineg only O- 

Therefore, each cell i is assigned the weighting coeffi- 
cient Aree Ws 
is just the cost of training a Marine in cell i normalized 
-l is the 


is the “training index" for cell i which 


by the maximum training cost over all cells. A; 


reciprocal of the actual number o£ Marines an cell 1 ee 
an exogenous weighting factor which permits the user of the 
model, if desired, to exercise his judgement in weighting 
cells or blocks of cells; presently, all We are SeE@Eo le 

The form of the resource function cD» is straighten 
ward. This function represents the expected number of 
dollars that must be allocated in order to offer cell i an 
SRB multiplier of x; and is the produck™@or the expected 
number of reenlistments in cell i that would result from 
effering abouts ae R,; (x, )E;; and the size of the bonus 
min ((./5)Po se, $12,000} where $12,000 is 75% £(the 
percentage of the bonus that the individual receives as a 
lump sum payment) of the maximum permissible bonus amount of 
5 16,0007 

Problem Pl is recognizable as a nonlinear knapsack 
problem in which cells represent commodities and the budget 
represents the weight or cube constraint. However, the 
nonlinearity in the objective function and budget constraint 


make solution by standard branch and bound techniques 


14 


impossible [Ref. 6], and make solution by dynamic program- 
ming inefficient Since Standard reduction techniques 
[Ref. 7] are inapplicable. A straightforward dynamic 
programming solution is possible using the cells as stages 
and dollars as the sole state variable. However, signifi- 
cant computational difficulties would arise due to the large 
number of stages (on the order of 10°) combined with a state 
variable that for FY86 can achieve any integer value on the 
interval [0,7x10°]. The solution methodology that is used 
could be applied directly to the nonlinear knapsack formula- 
tion. However, the explanation of the implementation of this 
methodology will be more transparent if the problem is first 
converted to a linear integer problem, in particular, a 
generalized assignment - problem. The following section 


adecertoesmenis mew Lormulation. 


B. CONVERSION TO A GENERALIZED ASSIGNMENT PROBLEM 
A typical formulation for the generalized assignment 
problem is: [Ref. 8] 


hmarces: 
ee as 1 agent 
eS VL | task 
Data: 
ee eest if task 1 is assigived ta. agent K 
eae amount of resource expended if task i is 
assigned to agent k 
CT ed minimum and maximum amounts of the 


resource that may be expended by agent k 
Decision Variables: 
eae 1 if task iis assigned to agent k 
Q otherwise 


Formulation: 


| 


iD 


yesip = gO re 1) 
k 


Xik € £0.15 


In contrast to the classical assignment problem, the 
first n constraints imply that more than one task i can be 
assigned for completion to agent k provided that the 
resource constraints on agent k are not violated. (Only one 
type of resource is allowed.) The second set of m 
constraints requires that each task be assigned exactly 
once. 

In the SRB multiplier selection problem the model is 
Simplified to the following: 


Indices: 
1 = lyon cells 
j = O, Sane SRB multiplier values 
Data: 
a eG 
os sy) 
B budget 
ne maximum permissible multiplier for cell i 
ge = 5,4,3 for Zone A,B,C, respectively 
Decision Variables: 
Xi j 1 if multiplier j is selected for cell i 
O otherwise 
Formulation: 


eg 
See Dat esas 8) 
ij 


Wes = 1) =e eee. 4m 


} 16 


sia ae a iO ab, 


Here, a Single agent is required to perform all m tasks but 
has at his disposal n; alternative processes for each. The 
above formulation is completely equivalent to Pl although a 
continuous, i.e., linear programming (LP), relaxation of 
this formulation might be very weak since neither Cx) nor 
re we typically ~be convex “functions. No problems 
arise, however, since no LP relaxation is employed in the 
solution procedure. 


This formulation of the model is quite general. The 


coefficients Cij and rij are directly calculated by evalu- 
ating the functions c,;(x;) and r;(x;) of the preceding 
FeocUloumacmealtesrcasible values of x, . Provided that the 


assumption of independence between cells is retained, alter- 
Mative functions are permissible. It will be seen in the 
succeeding chapter that the solution methodology is equally 
general and introduction of alternative functions is 


tc ie alee 
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TIL. SOLUTION METHODOLOGY 


With the problem and formulation now defined, this 
chapter details the approach taken to find a satisfactory 
set of SRB multipliers. Commonly, bounds for integer 
programs are established using an LP relaxation; the special 
structure of this problem and the likely weakness of the LP 
relaxation, however, suggests an alternative approach using 
Lagrangian relaxation. The relaxed formulation is extraor- 
dinarily easy to solve as an integer program and the 
resulting bound is at least as good as that produced using 
an LP relaxation. Furthermore, in the course of determining 
the bound, feasible sets of multipliers are encountered, the 
best of which is heuristically improved to yield a final 


Solution. 


A.  LAGRANGIAN RELAXATION TO ESTABLISH A LOWER BOUND 

The technique of Lagrangian relaxation is well suited 
for those integer programming problems that would be rela- 
tively easy to solve were it not for one or more compli- 
cating constraints. Moreover, it has been successfully and 
efficiently applied previously to generalized assignment 
problems [Refs. 8,9]. 

As described in Fisher's excellent primer on Lagrangian 
relaxation [Ref. 10], the approach is to move complicating 
constraints into the objective function using the product of 
the Lagrangian multiplier and the constraint violation as a 
penalty term. To demonstrate, consider the linear integer 


programming problem: 


pagh. (eek (P4) 
Shi COAX 9D 
Dx<e 


x>0O and integer 
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where the constraints Ax<b are the complicating ones. This 
is relaxed to: 

min, cx+)(Ax-b) (P5a) 

Sie. Dx <6 


x>0O and integer 


where it is assumed that this is easy to solve for a fixed 
row vector )A>0. Haeesstich a A, an optimal solution to 
problem P5 provides a lower bound on the objective value for 
P4. The best possible lower bound from a relaxation of the 


meminm Of PS 1S found by solving: 


max 5 9 min, cx+)(Ax-b) (P5b) 
Stax <e 


x>Q and integer 


Furthermore, this bound is at least as good as that provided 
by an LP relaxation [Ref. 11]. 

Examination of the generalized assignment formulation 
developed in the preceding chapter indicates that the only 
truly complicating constraint is the budget constraint; 


integer optimization with this constraint relaxed consists 


of m separable, "multiple choice" problems which are 
trivial. Thus, the relaxed formulation with scalar ) is: 
en 1 j 
Sc osieel 1 = lS aoe 
j 
ite Sally 


where the inner portion of the objective function may be 


equivalently written: 


min, yD 156i tarTiy)-AB (eqn a) 


Us) 


Ihe, 


For a fixed }, the inner minimization is performed by 
selecting a minimum (cj i FAQ 3) for each cell. The outer 
maximization problem iS a convex optimization problem which 
is easily solved, since } is a scalar variable, by first 
bracketing } and then performing a bisection search. 

The bracket about the optimal } is readily established. 
The lower bound iy)" is zero.) For sufi terently lass aa 
solution to the inner minimization is to expend as_ few 


dollars as possible, jase eee oc goes Ot 2 ae the upper 


u? 
bound on interesting values of ,, is the smallest value of 
for which this solution PS optimal. Therefore) (l= eee 


smallest value of } satisfying: 
CigtATI QS Ci 5 tATiS Vie (eqn 3.2) 


Recalling that Ti; represents the cost of offering cell i an 
SRB multiplier j and that j=0O corresponds to no bonus, Pope 


for all i. Thus, equation 3.2 becomes: 

Ci0SCj5tATiG V1. eee (eqn 3.3) | 
and, as a result, 

Ly=max{ (cyg-¢3 5) /r45} SM) Sa. (eqn 3.4) 


For )=L., 2 bie exists ber een the solution apm for valle 
and at least one solution where xij for some 1 and some 
j>1. For unusual data, the latter solution could be infeas- 
ible. To ensure that at least one feasible solution is 
obtained, in practice Ls replaced by MCL.) Peawhene lca sue 
The solution to P3b provides a lower bound on the 
optimal solution to the SRB multiplier problem. Upper bounds 
are provided by feasible sets of SRB multipliers encountered 


while solving P3b. Such solutions arise Since, for } 
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sufficiently large, the entire budget will not be consumed. 
The best of the feasible sets of SRB multipliers, that is, 
the one with the smallest objective function value and, 
correspondingly, the lowest upper bound value, is improved 
uSing the heuristic procedure described in the following 


section. 


be MARGINALS RATE OF SRETURN HEURISTIC 

It seems reasonable to expect that the best feasible 
solution obtained in the process of optimizing P3b would 
require the expenditure of nearly the entire SRB budget. 
There might, however, be some budget dollars remaining, a 
residual that could be apportioned to further reduce the 
objective value of the best feasible solution which is the 
present upper bound. Accordingly, a heuristic method for 
allocating this residual was developed using the concept of 
marginal rate of return. 

Given a best feasible set of SRB multipliers from the 
Lagrangian procedure, the incremental cost for each cell 
that would be realized by increasing that cell's multiplier 
by 1 is calculated for all those cells with multipliers not 
already at their maximum values. If this incremental amount 
is less than the unallocated portion of the budget it is 
possible to increase that cell's multiplier by 1 while 
remaining feasible. Let such a cell have x.;.;=l. Then, the 


13 
marginal rate of return for the cell is defined to be: 


ROR|=(€¢55-¢y  j4q)/ (84, j41°TGg) (eqn Jee) 


This/quotient expresses the improvement realized in the 
objective function per dollar spent when cell i's SRB multi- 
plier if increased. In the heuristic procedure, the cell 
with the maximum positive ROR; has its multiplier incre- 
mented’. The residual portion of the budget is reduced and 


the process is repeated until it is no longer possible to 


Zl 


increase the multiplier of any cell. The resulting set of 
SRB multipliers is accepted as the final solution and upper 
bound for the SRB multiplier problem provided that the ratio 
between the lower and upper bound is close to l. 

Other heuristic improvement procedures are possible but 
were not implemented in this study. For example, by simul- 
taneously incrementing one cell's SRB multiplier by 1 and 
decrementing another's (an O (m2) operation) 1t gmiemee be 
possible to improve on the upper bound while remaining 
feasible. In fact, the problem could be solved to optimality 
using a branch and bound algorithm [Ref. 9] though the large 
number of decision variables might make this a slow process. 
Pursuing strict optimality, with the resultant increase in 
algorithmic complexity and solution times was not deemed 
worthwhile in light of the consistently good solutions that 


have been obtained. 


ce 


IV. IMPLEMENTATION AND COMPUTATIONAL RESULTS 


A FORTRAN 7/7 computer program was written to implement 
the solution methodology described in the preceding chapter. 
Using data obtained from the Manpower Department at 
Headquarters, USMC (HQMC), the model was run on both an IBM 
3033AP and an IBM PC. To judge the robustness of the model, 
additional testing was performed uSing randomly perturbed 
data. 


A. COMPUTER PROGRAM 

The program that performs the Lagrangian bounding proce- 
dure and which yields the heuristically improved final solu- 
tion was written in ANSI standard FORTRAN 77. No functions 
unique to the system were used. When the program was run on 
the IBM 3033AP it was compiled by the IBM VS FORTRAN 
compiler at OPT(2); on the IBM PC it was compiled using 
Ryan-McFarland's RM/FORTRAN compiler with code optimization. 
The program requires six input files: Parameter file, 
Reenlistment Plan file, Reenlistment Rates file, Training 
Cost file, Cell Size file, and a Cell Weight file. File 
Specifications are detailed in Appendix D. As output, the 
program produces a file containing an MOS/Zone listing with 
the proposed SRB multipliers. For comparison purposes, the 
objective value as well as the percent deviation between the 
final solution's objective value and the lower bound on the 
optimal objective are output. 

Under some circumstances the uSer may wish to preassign 
Ehemorbemultipliers @in specific cells. For example, the 
uSer may wish to declare the multipliers for all three zones 
in MOS 2112 to be some specified value. This may be accom- 
plished by "tagging" those cells in the Cell Weight file and 


entering the desired SRB multiplier in the appropriate 


es: 


column in the data file. Additional details are provided in 
Appendix D. Comparison of the objective values before and 
after the introduction of these preassigned multipliers 
enables the uSer to gauge the effect of a przorz multiplier 


assignments. 


B. “GESi. Dale 

The FY86 data that was used to develop and test the 
model was obtained from the Enlisted Plans Section of the 
Manpower Department at HQMC. Model input consists of those 
values listed in equation 2.1 and 2.2 organized into five 
data files and one problem parameter file. Because the rele- 
vant data is extracted from a variety of sources, it was 
deemed desirable to place them in different files. The 
information contained within the Parameter file is primarily 
run specific. Samples from the various input files are 


contained in Appendix B. 


C. COMPUTATIONAL RESULTS 

The data described above leads to a problem with 979 
cells. This in turn results in a problem, in the form of 
P3b, with 4895 binary variables and 980 constraints sethe 
objective function value from the final solution is within 


.01% of the lower bound established through the Lagrangian 


procedure. In the process, 99.99% of the budget was 
expended. Appendix C contains a sample from the output 
file. 


Running time on the IBM 3033AP was 1.7 seconds. 1.4 of 
the 1.7 seconds is used to read in the data and calculate 
the coefficients Ci; and Tig: 
the heuristic use most of the remaining .3 seconds. The 


The Lagrangian procedure and 


time required to write the solution was negligible. 

After preliminary development and testing on the IBM 
3033AP, the program and data were copied onto a floppy disk- 
ette and, using the RM/FORTRAN software, installed on an IBM 
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PC configured with 512Kb RAM, an 8087 Math Coprocessor and 
two 360Kb floppy disk drives. No changes to the code were 
required in moving it to the PC. On the IBM PC, the program 
runs in approximately 4 minutes. Most of this 4 minutes is 
consumed reading in the data from the floppy diskette; 
Slightly faster times might be expected using a system with 
a hard disk arrangement. 

Additional testing was conducted to verify the robust- 
ness of the solution procedure and to confirm that the 
excellent results that were observed were not merely fortu- 
nate happenstance stemming from a "good" set of data. In 
this procedure, the data read by the model was perturbed in 


thLS manner: 


Geer UNITE (.5,1.5), 


where 
el) perturbed datum 
d original datum 
UNIF(a,b) random variate generated from the continuous 
uniform distribution between a and b 
LXJ largest integer not exceeding x 


This randomization process was performed on all the data 
except that response rates were capped at 1 and not discre- 
tized. 

In all, 50 additional model runs were conducted in this 
manner, each with a different data set. Solution times on 
the IBM 3033AP varied slightly, ranging between 1.5 and 1.9 
Seeconcdsemein each case, the final sol@Gion was within .027 
of the lower bound on the optimal solution. Typically, the 
heuristic improves the best feasible solution with respect 
to the lower bound by less than .04%. In a situation in 
which exceedingly tight bounds are not necessary, the 


heuristic could, in fact, be eliminated. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


In this thesis, the process of determining Selective 
Reenlistment Bonus multipliers subject to a limited budget 
was modeled as a mathematical program. In the course of 
doing this, an objective function based on expected devia- 
tions from cell reenlistment targets was defined. This is an 
important advance; formerly there was no method by which 
alternative sets of multipliers could be compared. Using 
this objective function, a procedure was developed which 
determines a nearly optimal set of SRB multipliers. The 
speed with which the solution is found, a modest core 
storage requirement and a compact, specialized FORTRAN code 
proves the procedure amenable for PC use. This isS an impor- 
tant advantage given the PC's wide-spread availability. The 
inclusion of an exogenous weighting factor and the capa- 
bility to preset multipliers allows the user additional 
flexibility in using the model. 

The model and solution method that are employed are very 
general and only require that the assumption of inter-cell 
independence be maintained. Other penalty and resource 
utilization functions are readily incorporated into the 
model. Improvements in the estimates for the response rates 
should further increase the accuracy with which the model's 
multipliers achieve the cell targets. 

Provided the user possesses the requisite hardware (the 
IBM PC configuration explained in Chapter IV should be 
considered the minimum) and a FORTRAN 77 compiler, the 
program is completely operational. A copy of the source 
code is listed in Appendix A. Other pertinent information, 
including examples of input and output files and user 


instructions, 1S contained within Appendices B, C, and D. 
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The future portends increaSingly restrictive budgets 
within the Department of Defense and additional pressure to 
maximize benefits given limited resources. With respect to 
maximizing the benefits realized from the SRB program, a 
nearly optimal assignment of SRB multipliers goes far 


towards achieving that goal. 
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APP 


ENDIX A 


LISTING OF SOURCE CODE 
IMPLICIT REAL* 8 H,O-Z 
REALS C(1100-0: “GA 0 
INTEGER*=2 MOS G5} .8 ZONE 858, IXBEST(1100) ,I1X(1100) 
COMMON CI CINE, BUDGET, pe NUMCE 
CALL INT TAL (c HLAM , MOS , IZONE, IX) 
GATEE Ie Bs ZU.1XBEST,HLAM,C,A,COST,MOS) 
CALL HEURI {Goa , LXBEST, ZH, “Gost. MOS)’ 
a Re ce 3 sess ves'e s'e Veveveslevevevevsveteletevevelele7%evc t ) 
ae site (1. D0 “OL ZH) 
FO Chk OS se CTIvE FUNCTION VALUE IS:',F11.3) 
WRITE (16 125 )PCTDEV 
‘OF (1X, OT LS SOLUMON el Se WiTHinws Stielec. 
/, HE OPTIMAL’ ) 
WRITE oe, 
DO 150 16 1 SRONCEL 
WRITE (16,200)MOS(L),IZONE(I1),IXBEST(L) 
FORMAT. (1X HS iad en als gO Rem) 
FORMAT tax Fa ek 3e> ZONE ,5X, MULTIPLIER’ ) 
STOP 
SUBROUTINE INITAL(C,A,HLAM,MOS,IZONE, IX) 
IMPLICIT INDECER®4 (12N)— 
REAL*8 AOS (3) sba¥(3) (1100 0:5),A(1100,0:5),R(0:5) 
INTEGER*2 MO (1) iWAe E 2 IX (1) 
OMMON CINF,BUDGET, TEST, NUMCEL 
READ (.10,100)BUDGET 
READ (10;200)NUMCEL 
READ (10;100 X 
READ (10;300)(AOS 1) rsil 3) 
READ (10,300)(APA ( ),1=1,3) 
READ (10; 100)BNSMA 
READ (10;400)PCTGFY 
READ (10,400)Q 
HLAM= -CINF 
CALE=1000.DO 
BUDGET=BUDGET / SCALE 
BNSMAX =BNSMAX*PCTGFY 
READ(11 73003 (05 (1) IZONE(T ik ACT, 0),I=1,NUMCEL) 
READ (14,6 00) (AC MO eAC MC 
ie pea dee: 
cone eae (159700) (C( ,J),J=0,MAXMUL ) 
READ (12,800)(A(I, 1B) 1 WL. oe 
READ (13,900)(A(1,4),1=1,NUMCEL 
DO 150 I=1,NUMCEL 
MAXMUL=6 - IZONE(TL) 
WI=A(T 0) 
E=A ( 1} 
TE ane } 
AGien (il een 
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WO 00 NIONUN Oo PO 
OOOO ODOOOOWWY 
ODODOOOO000W 


TI= sess 4) 

iaetl -.LT. 1.D0) Tr=50. D0 

DO Te een 

monn 82 

AVGPAY=APAY IZONE (I ) 

AVGAOS=AOS(IZONE 

IF (ACT .GT, 0.0D0)° THEN 
CW=WT* (TI /TIMAX) /ACT 


ELSE | 
CW=WT*TI /TIMAX 


END 
ieee ear. 1.D0) E=1.D0 


DO 125 J=0,.5 
I (J LE. Pe THEN 
IFTRM=TE-E=R 
IF (DIFIRM 100 DO) THEN 
~ eth Pe SCWEO (DI TRM**2 ) 


_ (1,3) =CH (DEFTRM2) 


BONUS= PCTGFYXAVGPAY*AVGAQS* J 
IF (BONUS GT. BN MAX) BO BONUS=BNSMAX 
ley )=R R(J)* . ‘BONUS /S AL 


al = (eT, 0) Cr AM}? eM oe ay 


C(1, 34 CINF 
A(I,J)=CINF 


AMBDA= ee ein 73 
CERO eae ahs ae on CUSE,©,4, MOS) 


ENDR=AMBDA 
153 vay ee a THEN 
UBEST=ZU 
FLAM=AMBDA 
END IF 
GOTO 100 
END IF 
ENDL=AM 
ENDR= 10. 0 1D0? “END 
AMBDA= (ENDL ENDR) /2.D0 
CALL MINFCN CAMBDA, ZL,1IX,ZU,COST,C,A,MOS) 
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100 


ine ne 
C CALCULA Se FEA 
CALL M 


100 


ae ore ae dy UII THEN 


ENDR=AM 
IF (20 aL ALE, PUB ESE) THEN 
FLAM= 2 AMBDA 
END IF 
GOTO 100 
Dir 
ENDL= ae 


DR 
AMB DA= (ENDL*ENDR [2250 
ZL, LX, ZU,COST ,C, Aymas) 


Cable of! Tre ABUB 
3 GET)’ THEN 
See 
IF (2 aL ALE, jZUBEST) THEN 
FLAM= S AMBDA 
END IF 


ELSE 
ENDL=AMBDA 
Deir 


NDR= eT . GI EPs) C070) 100 

IBLE SOLUTION 

FEA(FLAM, ZU, IXBEST,COST,C,A,1IX,MOS) 
ee 
SUBROUTINE ALS (AHO Zi eArRicOsl ,G,A,M0S) 
IMPLICIT REAL* on an 


IMPLICIT INTEGE N) 
INTEGERS? 1X i) “hoes (iMdo,0:1 
COMMON CINF, 1x(1) .MG her NUMCEL 
CELTOT= 6. DO 

ZU=0.D0 


DO 400 is 1 NUMCEL 
AMBDA*A(I,IX(1)) 


J=eCile J) + AMB Ea) 
CELOBJ2C sp 
IF (23 hE? ee ELMIN) THEN 
CELMIN=ZJ 
COBJ=CELOBJ 
INDEX=J 
END IF 
END IF 
CONTINUE 
CELTOT=CELTOT+CELMIN 
ZU=ZU+COBJ 
COST=COST+A(1, INDEX) 
CONTINUE 


ZL=CELTOT - AMBDA*BUDGET 
RETURN 


END 
SUBROUTINE MINFE 


A(FLAM, ZU, IXBEST,COST,C,A,1IX,MOS) 
IMPLICIT REAL*8 (A-H,0-Z) 
IMPLICIT INTEGER™4 (T-N) 
REAL* Csi 0: 1) A(LL 0,0:1 
RATEGERS IXBEST ( \,1X(1)3MOS(1) 
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COMMON CINF, BUDGET, TEST ,NUMCEL 
COST=0.DO 
CELTOT=0.DO 


ZU =0.D0 
DO 400 T=1 NYMC 
IF (MOS(I one ix? THEN 
Fee )) 
INDEX= 


X( 1 
pupa as aly) 
OTO 300 


100 hor200 J=0-5 
BEGG eeLE, TEST) THEN 
JECT: 1) FLAN “A(L,J) 


CELOBJ= J) 
EE Wauee LE. ELMIN) THEN 
LXBEST(T)=J 
COBJ=CELOBJ 
INDEX=J 
END IF 
END IF 
200 CONTINUE 
300 ZU=ZU+ 


COST= COSTA (I. INDEX ) 
400 CONTINUE 
RETURN 


END 
SUBROUTINE HEURIS(C,A, , ZH, COST ,MOS ) 
IMPLICIT REAL*8 (A- 
IMPLICIT INTEGER*4 
REAL*8 C(1100,0: 1) 
INTEGER*2 IXBEST 
COMMON CINF, BUDGET 
GOST L=COS 
100 RESID- BUDGET - cost l 


IF (MO te 0) THEN 


Bes} 
t-N) 
eet O0:1),UARRAY(1100) 
EST , NUMCEL 


I )=-CiINF 


(TXBEST(L) GE. 5 .OR. C(I, (IXBEST(I)+1)) 


* téq TEST) 4 
¥(1)=-CINF 
GOTO 
END IF 
RNUM=C(I are Gl ACL LAB ded 
RDENOM=A ne S 3"1dy B 
IF (RDENOM Aton sa 
UARRAY ( 
GOTO j? 
_YARRAY (I) =RNUM/RDENOM 
fe .LT. 0.DO) UARRAY(I)=-UARRAY (TI) 
200 CONTIN 
UBEST=--CINF 
DO 300 NUMC 


F (GA ARRAY ( L).GT. UBEST) THEN 
UBEST= “0 anh exh 


cary 


ELSE 


300 ea jiae Ne 


(UBE To TEST) ) GOTO 40 
GOs ie COSTA D EST { INDEX) ) 
*+A (INDEX, (IXB ST (INDEX) +1) 


SE 


400 


500 


IXBES ST (INDEX = [XBEST (INDEX 
IF (CO . BUDGET) GOT ‘Too 
ZH = 0.0DO 
COST=0.DO 
DO 500 I=1 
ZH=2i+C(T, 1 EST (I DI 
COST=CO IXBEST(I)) 
CONTINUE 
RETURN 
END 


a2 


APPENDIX B 
EXAMPLES OF INPUT FILES 


PARAMETER FILE 


70000000 
72000 : 
1010 29 1340 
16000 


in 


LO 
L300 
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REENLISTMENT PLAN FILE 
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CELL SEZE rime 
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TRAINING COST FILE 
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APPENDIX 


SAMPLE FROM OUTPUT FILE 
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APPENDIX D 
USER INSTRUCTIONS 


A. DATA FILES 

Input to the program is provided through six data files. 
File names,specifications and entries are detailed in the 
accompanying tables. In general, data values should be 
right justified within the fields. Rates, percentages, and 
the exogenous weighting factors may be placed anywhere 
within the specified field provided they are entered as two 
place decimals (e.g., .25, .03 or 1.15). With the exception 
of the Parameter file, all files contain the same entries in 
the first 8 columns. The program reads data values sequen- 
tially from the files. Therefore, it is vital that the 
MOSs, Zones and accompanying entries from one file corre- 
spond line for line with all other files. For example, if 
line 15 from one file contains entries pertaining to MOS 
0193, Zone C, then line 15 from all other files should refer 
to MOS 0193, Zone C. 

Those cells’ for which a non=ZeroueskS ) iu] tip iter ees 
expressly prohibited, such as the cells in the Band and 
Marine Corps Exchange MOSs, should be eliminated from the 
data files. Similarly, there are MOSs held only by very 
Senior Marines (E-8's and E-9's) not eligible for bonuses 
due to years-of-service requirements exceeding 14 years. 
These, too, should be purged from the data files. As 
written, SRBMULT.FOR can read a maximum of 1100 cells. If at 
some point it becomes necessary to increase this, it will 
require changing the array and matrix dimensions within the 


program from 1100 to the desired number. 
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Line 
& lal 


TABLE la 
DATA FILE SPECIFICATIONS 


Reenlistment Plan File (RPLAN.DAT) 


Column Data DeScription 


le» MOS 

b-7 blank 

8 Bomemc@letor Aw 2 fore, 3 for C) 
P= ES Number of Marines eligible for 


reenlistment in cell 
19-28 Reenlistment target for cell 


TABLE sib 
DATA FILE SPECIFICATIONS 


Cell Size File (ACTNUM.DAT) 


Column Data Description 


1-5 MOS 

6-7 blank 

8 mone Wluecorezam=. tor B, 3 for.C) 
9-14 Toca leanumber om Macines in cell 
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TABLES Le 
DATA FILE SPECIFICATIONS 


Trainingeeost Fale s( 1 eGsa na 


Column Data Description 


Zone (1 for 7A 2 for Bees fouee) 
Training cost’ (dollars} 


TABLE ld 
DATA FILE SPECIFICATIONS 


Response Rate File (RRATE. DAT) 


©) 
oO 
}— 
io 
= 
2) 


Data Description 


MOS 
blank 

Zone (1 for faeeeeoress. 3 for 
Response to SRB multiplier 
Response to multiplier 
Response to multiplier 
Response to multiplier 
Response to. SRB multiplier 

if sep hiee ote 

esponse to SK 
(if applicable) 


| 
in 


WW Wher wWodOd+ 
Oo WNW! 
Piet bees 
WWhDYOL 
COMA 
— 


i 
-—* 
J: 
Wm PwWMROXM 


multiplier 
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TABLE le 


DAEA PrebeoeeECiIFICATIONs 


Cell Weight File (WEIGHT.DAT) 


Data Description 


MOS 

blank 

Zonecmi@letOr A. 2estomeb., 39 tor C) 
lank 

Preset uokbemultrp lier (Ovilt not 
Bese 
lank 

Exogenous Soe oan factor (typic- 

ally, use 1.00. If other values are 

used, include decimal, e.g. 0) 


TABLE lf 


DATA FILE SPECIFICATIONS 


Column 


‘a 

Crh ot CN ea 

PRN FP RENE oe 

mM ME NE MH NHN 
iA Re 


bE RE COR 100-4 


- Fr 


Parameter File (PARAM.DAT) 


Data Description 


Budget available for current year 

bonuses (dollars 

Total number of cells 

Maximum EE eLE RTE GOR: (dollars) 

Zone A avg reenlistment duration 

Zone B avg reenlistment duration 

Zone C avg reenlistment duration 

years 

Zone A avg monthly base pay 

Zone B avg monthly base pay 

Zone C avg monthly base pay 
ectaLilaiee 


Maximum permissible bonus (Ce ERS) 


Percentage of bonus awarded as_l 
sum en ecimal form, e.g., ./5) 
Over/under factor 


ump 
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B. RUNNING THE PROGRAM 

To use the program, the user should possess, at a 
minimum, the IBM PC configuration detailed in Chapter IV and 
a FORTRAN 7/7 compiler. The program was developed and tested 
uSing the Ryan-McFarland RM/FORTRAN compiler. However, the 
code used in the program is sufficiently generic that other 
FORTRAN 77 compilers should be acceptable. 

A diskette ("program diskette") containing the following 
files will be provided to the Enlisted Plans Section, 
Manpower Department, HOMC: HOWTO.DOC, RUNCHK. BAT, 
FILGCHKSEOR; PHeGHK EXE. RUNSRB. BAT, SREMULT . FPORe and 
SRBMULT. EXE. The file HOWTO.DOC contains a copy of oper- 
ating instructions. Files with filetype .BAT are executive 
routines, those with filetype .FOR are FORTRAN source code 
and files with filetype .EXE are executable programs. Use of 
these programs is detailed in the following Sections. Prior 
to running any of these programs certain parameters within 
the DOS configuration file CONFIG.SYS must be changed. If 
this file does not already exist, it will be necessary to 
create it. The RM/FORTRAN that was used to compile 
SRBMULT.EXE and FILCHK.EXE requires that the number of files 
that can be opened concurrently be increased to at least 
ten. Therefore, the following line must be added to 
CONFIG.SYS: 


IE as Ly, 


Similarly, the number of disk buffers allocated by DOS at 


system startup must be increased to ten. Accordingly, add: 
BUFFERS = 10 


to CONFIG.SYS. Once these changes have been made_ to 
CONFEG . SYoumrebooreDes. 
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i[eecnecorne Datanri les tor Conformity 
The FORTRAN program FILCHK.EXE is provided to check 


Bice umes ber MOS and Zone conformity prior to using 
SReEMULT. FOR. To use FILCHK.EXE on a system with two disk 
drives, place the program diskette in B drive and the disk- 
ette containing the data files ("data diskette") in A drive. 


Enter the following DOS commands: 


Pace = B2\ 
A: 


The first command will cause DOS to check the B drive for 
executable files. The second sets the A drive as the default 
drive. At the system prompt, enter RUNCHK to execute the 
program RUNCHK.BAT. This executive routine will Set neces- 
sary file definitions and then execute FILCHK.EXE. On a 
system with a hard disk, place all files in the same direc- 
tory on the C drive (hard disk) and then execute RUNCHK.BAT 
as before. 

The FILCHK program uses the Sequence of MOSs_ and 
Zones in the file WEIGHT.DAT as a reference set and checks 
Ehewotmermrour data files for conformance to it. If a non- 
conforming entry is discovered, the applicable file defini- 
tion number and line number are written to the Screen. The 
remaining files are then checked. Once detected errors are 
corrected, RUNCHK.BAT should be executed again until all 
files conform with the reference aS evidenced by no error 
Perirns: 

2. Running SRBMULT 

To run SRBMULT.EXE on a system with two disk drives, 

place the program diskette in the B drive and the data disk- 


ette in the A drive. Enter the following DOS commands: 


Pa ihe 5° \ 
A: 
RUNSRB 
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Several minutes later, an output file MULPLN.DAT will be 
created on the diskette in the A drive. It will contain an 
MOS and Zone listing, the recommended SRB multipliers, the 
objective function value and the percent deviation from the 
optimal. The objective function value is useful in judging 
the effect of introducing exogenous weighting factors or 
presetting multipliers. Instructions on how to do both of 
these is contained in the following section. 

Running the program is somewhat simpler on a system 
with a hard disk. All files should be placed in a single 
directory on the C drive (hard disk) and that drive should 
be made the default drive. SREMULT.EXE “is inievabved ee by 
typing RUNSRB to execute RUNSRB.BAT. 

C. EXOGENOUS WEIGHTING FACTORS, PRESET MULTIPLIERS AND THE 

OVER/UNDER FACTOR 

Inclusion of an exogenous weighting factor (columns 
14-18 in WEIGHT.DAT) permits the user to exercise additional 
control in weighting the effects of reenlistment shortfalls 
in particular cells based on his judgement. Weighting 
factorswWi greater than the default alucsot .) 0 vesmheaeen 
comparatively greater penalties for failing to meet reen- 
listment targets. Conversely, weighting factors less than 
1.0 result in lesser penalties. The effect of the factor is 
linear. Thus a 2.0 results in a penalty twice as large as a 
1.0 and a factor of .5 reduces the penalty by one-half. As 
a weighting scheme based on training costs and manning 
levels is already utilized, this capabililty should be used 
with caution. In general, a value of 1.0 should be used for 
alli@eelis:. 

The user can also preset SRB multipliers. This might be 
necessary if factors other than thoSe considered by the 
model dictate that particular cells be assigned specific 
multipliers. This is accomplished by entering a minus sign 


(-) in column 1 of the pertinent cell's data entry in the 
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file WEIGHT.DAT. The desired multiplier is then entered in 
column ll of the same line. In the output file, the minus 
Sign will appear next to the appropriate MOS/Zone thereby 
“flagging” those cells with predetermined multipliers. 

The last entry in PARAM.DAT is the over/under factor Q, 
the weighting factor employed to permit the user to estab- 
lish the relative importance of overshooting or under- 
Shooting reenlistment targets. For example, a over/under 
factor of .70 implies that the penalty incurred in the model 
for overshooting the reenlistment target by some number is 
only 704 of the penalty that would be incurred for under- 
Shooting the target by the same number. Since an overage in 
a particular cell is generally considered more acceptable 


than a shortage, the factor will typically be less than l. 
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